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Abstract 

In this paper we study a phase structure of 5 D AT = 1 super Yang-Mills theory with 
massive matter multiplets and SU(N ) gauge group. In particular, we are interested in 
two cases: theory with Nf massive hypermultiplets in the fundamental representation 
and theory with one adjoint massive hypermultiplet. If these theories are considered 
on S 5 their partition functions can be localized to matrix integrals, which can be 
approximated by their values at saddle points in the large-A limit. We solve saddle 
point equations corresponding to the decompactification limit of both theories. We 
find that in the case of the fundamental hypermultiplets theory experiences third-order 
phase transition when coupling is varied. We also show that in the case of one adjoint 
hypermultiplet theory experiences infinite chain of third-order phase transitions, while 
interpolating between weak and strong coupling regimes. 
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1 Introduction and Main Results 

Recently, 5 D super Yang-Mills (SYM) theory attracted much attention due to its relation 
with 6 D (2,0) superconformal theory discussed in (l|. Using localization [2| it is possible to 
reduce full path integral of 5 D SU(N ) SYM to a finite-dimensional matrix integral |3,4|. The 
last one appears to be solvable in certain limits. In particular the free energy of 5 D JV = 1 
SYM with the adjoint hypermultiplet was derived using localization in [5]-[7|. It was shown 
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that the free energy of this theory behaves as N 3 in the strong coupling limit. This result is 
consistent with the well-known N 3 behavior of the free energy of 6 D theory obtained from 
the supergravity considerations |8|. Localization results were later generalized to theories 
living on the manifolds with more complicated geometries in [9-11 , revealing the same N 3 
behavior of the free energy with only difference in the prefactor. 

Another 5 D theory of interest is the super Yang-Mills with the USp(N) gauge group and 
infinite coupling constant, which is conformal fixed point in five dimensions. This theory is 
especially interesting because it has known holographic dual in AdS§ space |l2|. This gauge 
theory was studied using localization as well in 13 and 14 . It was found that the free 


energy of this SYM theory has N 5//2 behavior in the planar limit. This result matches the 
result obtained from the AdS/CFT correspondence. Finally this 5 D CFT was also studied 


on the backgrounds with more general geometries, e.g. on the squashed spheres, in 15[ 16 
These results allow one to evaluate supersymmetric Renyi in five dimensions 
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One can also introduce a Chern-Simons term into SU(N ) SYM theory and consider the 
limit of the infinite Yang-Mills coupling. Then this theory is also 5 D superconformal fixed 
point, which was studied using localization in 18 . It was found that in this case the free 


energy of theory also behaves as N 5 / 2 when N is large. This suggests the existence of the 
holographic dual of 5 D supersymmetric Chern-Simons theory, though the dual is not known 
yet. 

Using localization, it was recently found that different supersymmetric theories have 
a very interesting phase structure. The matrix models obtained from the localization of 
supersymmetric hied theories with the massive hypermultiplets were shown to experience a 
phase transition at some critical values of the couplings in the decompactihcation limit. First 
evidences of these phase transitions were obtained for 4 D J\f = 2* SYM on S' 4 . It was found 
that in the limit of infinite radius this theory experiences infinite chain of phase transitions 
as coupling is varied from weak to strong. In particular, the authors of 19 solved the matrix 
model obtained from the localization of Af = 2* SYM. They have found that the support of 
the eigenvalue density in this case consists of many intervals of the length equal to the mass 
m of the adjoint hypermultiplet. On the boundaries of these intervals density distribution 
has cusps. When the ’t Hooft coupling is increased the length of the support grows as well 
and at some points new cusps appears in the distribution. Emergence of new cusps signals 
about the phase transition at the corresponding point. Finally, at very strong coupling cusps 
of the distribution smooth out and the solution approaches the one obtained in 20 . This 
interesting phase structure of Af = 2* SYM on S' 4 was investigated in details in the series 
of papers 19,21,22 . Later these results were generalized to the decompactihcation limit of 
Af = 2* SYM on the ellipsoids in 123 . 

Another AD theory experiencing large -N phase transition is Af = 2 SYM with 2 N mas¬ 
sive hypermultiplets in the fundamental representation considered on S' 4 
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However, 

phase structure of this theory appears to be much simpler. It experiences only one phase 
transition, while interpolating between weak and strong coupling in the decompactihcation 
limit. From the matrix model point of view this phase transition takes place when the length 
of the support of the eigenvalue distribution becomes larger than two times the mass of the 
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hypermultiplet. 

Similar phase structures were also obtained in 3D supersymmetric theories. In 
authors considered matrix model obtained from the localization of Af = 2 supersymmetric 
Chern-Simons theory on S 3 coupled to 2Nf massive hypermultiplets in the fundamental 
representation. It was shown that in the decompactihcation limit this theory experiences 
third-order phase transition similar to the one in Af = 2 SYM with 2 N massive fundamental 
hypermultiplets on S' 4 . Finally these results were generalized to the case of massive defor¬ 
mations of the ABJM theory on S 3 26,27 , for which the phase structure of the theory 


appears to be similar to the case of Af = 2* SYM on S' 4 . Namely, it was found that the 
theory experiences an infinite chain of third-order phase transitions, while interpolating be¬ 
tween weak and strong coupling regimes. As in 4 D case, these phase transitions are related 
to the emergence of new cusps in the distribution of the eigenvalue density. 

In this paper we try to generalize the results of 3D and 4 D to the case of 5 D gauge 
theorie^] We are particularly interested in the SU(N) Af = 1 SYM coupled to either one 
adjoint or Nf fundamental massive hypermultiplets. To work with these theories we use 
the matrix model obtained by localizing AT = 1 SYM on S 5 with arbitrary hypermultiplet 
content |3j,j4] . In general it is not possible to evaluate these matrix integrals. However in this 
paper we focus only on the large -N limit where the matrix integrals are dominated by their 
saddle points. In the case of adjoint hypermultiplet this saddle point equations have been 
solved in the weak and strong coupling limits |6,7|. Here we send the radius of the five-sphere 
to infinity, which will simplify the saddle point equations and allow us to solve them exactly 
at any coupling. When the solution is observed we can evaluate supersymmetric observables 
in AT = 1 SYM. In particular, we consider the free energy and the expectation value of the 
circular Wilson loop. Studying both these observables we can draw conclusions about details 
of the phase structure of the theory. 

This analysis leads to the following picture of the phase structure of considered theories. 
In the case of the theory coupled to Nf fundamental hypermultiplets there is only one third- 
order phase transition taking place at 


t, = 


8n 


Nf 

1 - L 

m \ 2 N 


-i 


( 1 . 1 ) 


where t = g\ M -A is the definition of the’t Hooft coupling used throughout in this paper and 
m is the mass of the hypermultiplets. Notice that in order to reach this critical point one 
should consider negative g\ M - As explained in section [2] this does not contradict anything. 
Negative coupling arises in the theory due to the renormalization and doesn’t spoil the 
convergence of the matrix integral. 

The phase structure of the SU(N) AT = 1 SYM with the adjoint hypermultiplet is even 
more interesting. For this theory we observe an infinite chain of third-order phase transitions 


1 Some phase transitions were already found before in 5 D J\f = 1 Chern-Simons theory |l8j and in 5 D 

However these phase transitions are different from the 


Af = 1 SYM theory with adjoint hypermultiplet 28 
transitions described in this work. 
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at the following critical points 


Rtt 2 

*<"> = — (n+1) 

m 


( 1 . 2 ) 


In the strong coupling limit these phase transitions smooth out and approach the strong 
coupling solution found in |7|. A nature of these phase transitions is the same as in the 
mass-deformed ABJM theory and in 4 D AT = 2* SYM. As we show, the eigenvalue density 
at the saddle point of the matrix model has cusps. The number of these cusps grows as we 
increase the ’t Hooft coupling, and each time a new cusps appear the system experiences a 
phase transition. 

This paper is organized as follows. In section 2 we review some properties of the matrix 
model obtained by localizing 5 D SYM on S 5 . In particular we consider issue of renormal¬ 
ization of the coupling and find the decompactihcation limit of this model. In section 3 we 
solve the saddle point equations of the matrix model corresponding to the theory coupled 
to Nf fundamental hypermultiplets in the decompactihcation limit. Then calculating the 
free energy of the matrix model we show that the theory experiences the third-order phase 
transition at the critical value of the coupling. In section 4 we solve the matrix model with 
adjoint massive hypermultiplet in the decompactihcation limit and show how the chain of 
phase transitions described above arise in this case. Finally, in section 5 we discuss the effect 
of the squashing of hve-sphere on the phase transition. We find that in contrast to the AD 
case discussed in 23 , phase transitions are not affected by the deformations of the sphere. 


Some technical details of calculations are included in the appendicies of the paper. 


2 Matrix Model 

To study the phase structure of 5 d SYM we use the result of supersymmetric localization [3, ! 4], 
that reduces full held theory path integral to the hnite-dimensional matrix integral given by 




9 ym 


-Tr(^ 2 )-YTr(<^ 3 ) 


+ 0(e 


( 2 . 1 ) 


Here gyM is Yang-Mills coupling, k is Chern-Simons (CS) level, r is the radius of S 5 and 
integration variable </> is related to the expectation value of the scalar held a from vector mul- 
tiplet as the following </> = — zraj^J One-loop contributions of the vector- and hypermultiplets 
are given by 


zScpO= nn h - y 

P 0 


(l+ft+lt 2 ) 


( 2 . 2 ) 


2 Note that in order to have well defined matrix integral we should integrate over real <f> or equivalently 

imaginary a. In all our conventions we follow |4j. 
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and 


z i- p ioo P (4>) = nn ( t_ ^^ _m+ 2 


-(l+§t+H 2 ) 


(2.3) 


where /3 are the roots, /i are the weights of the representation R and m = —iMr with M 
being the mass of the the hypcrmltiplet. The matrix model described above was well studied 
in some regimes. In particular, planar limit of the SU(N ) SYM was studied in (6,7,29 . 
In 13,(14 this model was studied for the case of 5 d superconformal theory. And, finally, 


planar limit of the CS theory was studied in 18 


In general (2.1) should also include instanton contributions. However we, in this paper 


will be interested in the large-IV limit of the theory, which suppresses instanton contributions. 
Hence we omit all instanton contributions and consider only classical action together with 
one-loop determinants in the partition function. 

In this paper we put CS term to zero and concentrate on the decompactihcation limit 
of the SU(N ) SYM with different hypermultiplet content. In all cases hypermultiplets are 
massive which leads to the phase transitions similar to the ones obtained in [T9, 26, 2l]. 
Before moving to the solutions of the matrix model we should discuss general properties of 
this model. More detailed analysis of this matrix model can be found in |7|. 


2.1 Renormalization of the Coupling Constant 

As we can see the one-loop contributions for both vector- and hypermultiplets can be written 
in the following form 


zrLpM = -Y 

(2.4) 

z \ y -iZM) = V 1 ((/b , 

(2.5) 

where V(x) is the infinite product 


OO 

V(x) = x JJ (t - x) {1+ ^ t2) (t + x) (1 -§ t+ ^ 2) . 

t=l 

(2.6) 

This product is divergent with the divergent part given by 


°° ( x 2 \ 

log V = ( -3a; —— j + convergent part, 

(2.7) 


where we have also omitted x -independent divergent part, 
we introduce the UV cut-off at t 0 = ttA 0 r leading to 


To regularize the product (2.6) 


logV 



( 2 . 8 ) 
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Using this expression we extract divergent terms in the one-loop determinants 


log zr_i 0 M) 


2((/M0» 2 + 3 (/3u0) 


-vrA 0 r Y 

0 L 

7tA 0 r C , 2 (adj)Tr(^ 2 ) + convergent part 

lo g^o e o P (0) = 71A 0 rY 


+ convergent part = 


(2.9) 




+ 


convergent part = —tcAqt C*2(i?)Tr(0 2 ) + convergent part (2.10) 

Here C- 2 (A) is the quadratic Casimir element for representation R , dehned by Tr(T^T b) = 

C 2 (R)Sab- With these conventions XX (‘/bh)) 2 = 2C*2(i?)Tr(0 2 ). Note that we have omitted 

v 

terms linear in 0, because the gauge group is SU(N). In this case the tracelessness condition 
implied for (j) forces all linear terms to be zero. 

The divergent terms in the one-loo p co ntributions (2.9) and (2.10) are proportional to 
Tr(0 2 ). So in the full matrix integral (2.1) these divergent terms can be absorbed into the 
renormalization of the Yang-Mills coupling 


9eff 


9ym 


~ & I 


dj) - , 


( 2 , 11 ) 


where the sum in the second term is over all hypermultiplets of the theory. Equation (2.11|) 
reproduces the renormalization obtained in flat space using perturbation theory in 30 


Further we will consider theories with particular matter contents, leading to different 
effective coupling constants. In each case we will come back to (2.11) and consider particular 
examples of this renormalization. 

After the regularization described above finite parts of the one-loop contributions (2.9) 
and (2.10) can be written in the form 


log = -tiAi F v (<t>) , 

iogd'C” 8 ) w^-E tr «. 


( 2 . 12 ) 

(2.13) 


where Fv(4>) and F H {(f)) are functions dehned as follows 
F V (x) : 


“ log (sinh 2 (7nr)) - ^f(ix), 


1 


F h (x) = -li¬ 


nn 


ix ) H —l ( —b im + ix 


+ - f [ - — im — ix ) + - f [ - + im + ix 


(2.14) 


(2.15) 
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Here functions f{x) and l(x) are the ones introduced in |3| and 31 correspondingly. For the 


definitions and properties of these functions see appendix [A} Using results of the regulariza¬ 
tion it is convenient to write down the matrix integral ( 2 . 1 ) in the following form 


Z = 


[d(p]e 


-n® 


Cart an 

87r 3 * r 


= — 5 —tr + tr Ad F v (4>) + J^tr Ri F h (4>) ■ 

g eff v J , 


(2.16) 


Another observable that can be evaluated using localization technique is the expectation 
value of the circular Wilson loops 


(W(C)) = ( -Tr V exp 


dr (■ iA + ojx|) 


(2.17) 


where C is wrapping the equator of S 5 . Such loops were considered in 1321 and |7| for 
5 D SU(N) SYM theory and in [14] for 5 D superconformal theories. After the localization 


expectation value (2.17) in the case of fundamental representation turns into the expectation 


value of the matrix model (2.16) (2 


{W) = -(Tr <?**). 


(2.18) 


2.2 Decompactification limit 


The matrix integral (2.16) is not solvable even in the planar limit. Some particular limits 

for the case of one adjoint hypermultiplet and Nf 


of this matrix model were studied in 
fundamental hypermultiplets. In this paper we study the behavior of the theory with SU(N) 
gauge group and different hypermultiplets content in the decompactification (infinite radius) 
limit. 

I 11 all equations above we have used the dimensionless variables 0 and m, or equivalently 
we can assume that the radius of the sphere was set to one. To reproduce the r-dependence 
of the matrix integral we should rescale our variables as follow^} 

0 —> (f>r , fh —)■ mr . (2.19) 

The decompactification limit then corresponds to the limit r —> 00 . Using asymptotic 


expressions (A. 6 ) we can find that in this limit F(0) in (2.16) can be written as 

1 X^’( 0 )= ^ 

g eff 


2irr 3 


47T 2 / 1 1 \ 

F (0) = 0 2 ) + tV Ad - 2^2 101J 

-J2 trR i (^l0 + m l 3 + ^l0 + m l) 


( 2 . 20 ) 


3 From now on we will denote all dimensionless parameters with tilde and corresponding dimensionfull 

ones - without tilde 















Here we also included 1/r 2 terms linear in 0. These terms will be important when we discuss 
the effects of finite r in the large radius limit. 


We see that in the decompactihcation limit the matrix model simplifies a lot, as (2.20) 


does not contain any complicated functions and, as we show further, can be solved explicitly 
in the planar limit. However, this is not the only reason to consider this limit. As mentioned 
in the introduction, studies of the decompactihcation limit in different 3 D and 4 D theories 
revealed interesting phase structure of these theories. For instance, phase transitions were 
found in 3 D Chern-Simons coupled to the massive fundamental hypermultiplct 24 , in AD 
Af = 2 super-QCD with massive quarks in the Veneziano limit 21 , in mass-deformed ABJM 


theory 26,27 and finally in AD Af = 2* SYM theory [19], In each case theory was considered 


in the decompactihcation limit and it was shown that in the finite volume phase transitions 
disappear. We expect to observe similar picture in 5 D SYM theory and hence we should 


also concentrate on the decompactihcation limit of the corresponding matrix integral (2.16). 

All considerations above are general, 
hypermultiplet content 


They can be applied to any gauge group and 
From now on we will concentrate on the SU(N) gauge group, 
meaning that C 2 (adj) = N and C 2 (fund) = Inspired by the previous results on the phase 
structure of the supersymmetric gauge theories we consider two type of theories: Af — 1 
SYM coupled to Nf massive fundamental hypermultiplcts and Af = 1 SYM coupled to one 
massive adjoint hypermultiplet. 


3 Af = 1 SYM with Nf fundamental hypermultiplets 

In this section we consider Af = 1 SYM containing vector multiplet and Nf hypermultiplets 
with the mass m. Half of the hypermultiplets is in fundamental, while another half is in 
antifundamental representation!]] Before working out matrix model equations of motion 
and its solutions we should discuss the renormalization issues, as they play very important 
role in this case. According to (2.11) effective coupling in the case of Nf fundamental 
hypermultiplets have the form 


9eff 


9o 


Ap N 

8 i r 2 


¥ 


(3.1) 


where ( = Nf/N is the Veneziano parameter. Note that the theory has constraint on f 
that can be found from the decompactihcation limit expression (2.20). We want F(0) to be 


positive at large 0 in oder for the matrix integral (2.16) to be convergent. In the case of the 


SU(N) gauge group with Nf/2 fundamental and Nf/2 antifundamental hypermultiplets we 
obtain: 


2nr 3 


F (</>) = (l^ + m l 3 +1& _m l 3 ) + '' ' = 

j¥=i ' 




N,) J2 i* 


(3.2) 


+ • • • 


4 For shortness we often address this case as the case of Nf fundamental hypermultiplets in this paper. 
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where by the ellipses we mean all terms subleading in 1/|0|. From this expression we see 
that in order to have well defined matrix integral, parameters of the theory should satisfy 
Nf < 2N or ( < 2 in terms of the Veneziano parameter. The same restriction was obtained 
from the flat space prepotantial in [33]. 

According to (3.1) YM coupling can be negative, except in the case of £ = 2. Further 
in our considerations we will always assume ( < 2 . Then we can redefine UV cut-off A 0 so 
that 


9eff 


AN 

8n 2 


(3.3) 


and use it in the matrix model further. 

Now we are ready to solve the matrix model with described matter content and coupling 
renormalization. In the planar limit the matrix integral (2.16) is dominated by the saddle 
point satisfying 


dF((p) d ( 8 tt^t v—' ~2 v—v v—v „ ,7 

vr = 3i w&HEFw- 

dcj) k dct) k \ g eff V & i 


,) + ^ Y ( F «U) + 


= 0 . 


(3.4) 


Using expressions (|A.5|) for the derivatives of F v and F H we obtain the saddle point equations 

T 

+ (d>; + mY 

167r 2 r 


^ ^2 - (4>i - fa ) 2 j coth(7r(0i - fa)) + ^ Q + (fa + m) 2 j tanh(7r(0j + m))+ 

■jN V ' ' 

r /1 \ ^ Wn 2 r ~ ~ 

- ( - + (<t>i - fn) 2 ) tanh(7r(0j - fnj) = — & = -2A Nr fa . 

: V 4 / g eff 


N 

~4 


(3.5) 


After taking the decompactihcation limit we arrive to 

- 2 A Ncfii = ^Nf(<f>i + m) 2 sign( 0 j + m) + - m) 2 sign( 0 j - m) — 

- 0 j) 2 sign( 0 * - #,-). 

This equation can be also observed from ( |2.20[ ). Introducing usual eigenvalue density 


(3.6) 


(3.7) 


we rewrite the saddle point equation (3.6) as the integral equation 

—2A cf> = ~C(0 + m) 2 sign(0 + m) + ~C(0 * ?n) 2 sign((/> — m) — 

/ <%#)(0- V’) 2 sign(0 --0). 


(3.8) 
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We assume that the eigenvalue density p(0) has finite support [—a, a}. Then we should 
consider two separate cases. In the first case, which corresponds to the phase of the theory 
below the transition, a < m is satisfied. In the second case we have a > m, which corresponds 
to the phase above the transition point. 


To solve equation (3.8) we differentiate it three times: 


— 2A = —C(0 + m)sign(0 + m) + -</(0 — m)sign(0 — m) — 


1 

2 ( 


1 

2 ( 


2 J d0p(0)(0-0)sign(0-0), 

(3.9) 

b — m) — 2 J d0p(0)sign(0 — -0), 

(3.10) 

C (5(0 + m) + 5(0 - m)) - 4p(0), 

(3.11) 


where (3.9),(3.10) and (3.11) correspond to the first, second and third derivatives of the 


equation (3.8) respectively. Now we consider these equations for the cases a < m and a > m 


separately. 

3.1 Solution below the transition point 


When the theory is below the transition point (Phase I) terms with masses in (3.10) cancel 
each other and we are left with 


0 = / dif) p(-0)sign(0 — -0) = / d0p(0) — / d0p(0 ). 


(3.12) 


This equation should be satisfied for any 0 e [—a, a], hence we conclude that p(0) is sharply 


peaked near the endpoints of the distribution at 0 = ±a. Then the solution to (3.12) is 
given by 


= t: (W - a) + <5(<A + a)) , 


(3.13) 


where the overall coefficient 1/2 comes from the standard normalization condition 


d0p(0) = 1. 


(3.14) 


To be more precise the arguments of the 5-functions in (3.13) should look like (0 ± a =ft) 
with e —> 0. This small shift should be done in order to include the whole delta function 
into the density support [—a, a]. Everywhere further we imply these shifts of the 5-function 
arguments at the distribution endpoints. 


Note that (3.13) does not satisfy (3.11) at the endpoints of distribution. The best way to 


make the solution (3.13) more rigorous is to consider equations of motion (3.5) in the limit 
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of large but finite r, taking into account terms subleading in 1 jr. These terms correspond 
to the repulsion of the eigenvalues at small separations, which washes out (5-functions and 
turn them into the peaks of 1/r width. The form of these peaks can be found analytically. 


Then the solution (3.13) can be obtained in the limit r —y oo. Details of these calculations 
can be found in appendix [Cj 


To determine the position of the endpoint a, we substitute solution (3.13) back into the 


original equation (3.8) which leads to 


a = A + -(m 
2 


(3.15) 


By construction the solution (3.13) works when a < m or A < m (l — Notice that the 
r.h.s of (3.15) is always positive due to the condition ( < 2, we have discussed previously. 


3.2 Solution above the transition point 

Now we increase A, pass the point A c = m (l — and arrive to the phase above the 
transition point (Phase If), where we have a > m. To solve (3.8) we again use its three 


derivatives (3.9)-(3.11). First we address (3.11). In the case a < m terms S((j) ± m) do 


not contribute into this equation, because arguments of these h-functions are never zero. 
However, if a > m h-functions start contributing into the eigenvalue density. Then it is 


natural to assume that solution contains these h-functions on top of the solution (3.13) 


P iH) W = - [(2 - CW + «) + W + m) + W - m) + (2 - <)<# 


a)] . (3.16) 


Here the coefficients in front of <5(0 ± m) are chosen to satisfy (3.11) and the coefficients in 
front of <5(0 ± a) are found from the normalization condition (3.14). 


Now we substitute the ansatz (3.16) into original equation of motion (3.8) in order to 


obtain the position of the distribution endpoint a 




-l 


A. 


(3.17) 


This solution works when A > m (l — . Thus we covered all the values of A and, unlike 

in 3D Chern-Simons with Nf fundamental hypermultiplets 24 , we do not observe any 


intermediate phase. The only transition takes place at the critical point 


A c = m 




(3.18) 


The solution found here suffers the same problems as (3.13). To justify it we again consider 


equations of motion with large but finite r and take the decompactikcation limit in the very 
end. This calculation is done in appendix [Cj where we reproduce the eigenvalue density 


(3.16) in the limit r —> oo. 
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P(0) P(0) 



Figure 1: Eigenvalue density p(0) calculated with the following value of the parameters: 
r = 10, m — 5, C = N — 100 corresponding to A c ps 3.76. The orange dots show the 
results for the numerical solution, while the dashed blue lines show the positions of the 
h-functions in the analytical solutions (3.13) and (3.16). 


We also check solutions (3.13) and (3.16) numerically. In FigjTlthe orange dots show the 
results for the numerical solution of the full equation of motions ( |3.5[ ). On the same plots 
the dashed blue lines denotes the positions of the (5-functions in our analytical solutions for 


the densities below (3.13) and above (3.16) the transition point. As we see from these plots, 


our analytical solutions reproduce the result of the numerical evaluation very well. For some 
technical details of the numerical simulation see appendix [B] 


3.3 Free energy and the order of transition 

We now determine if the transition between phases I and II is phase transition and , if so, 
what is the order of this transition. To answer these questions we find the behavior of the 
free energy near the critical point A c . To calculate the free energy in the decompactification 
limit we directly use the asymptotic formula (2.20) with an appropriate matter content. In 


the large -N limit we express the free energy in the integral form 


27rr 3 A r2 


F = 


A J d0p(0)0 2 + 
c 



d0d#(0)p(VOI0~ Vf 


(3.19) 


— / d0p(0) (|0 + m\ d + 10 


m\ 


These integrals can be evaluated both below and above the phase transition point using cor¬ 


responding expressions (3.13) and (3.16) for the eigenvalue density. After some calculations 
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we obtain 


pi 1 ) — 
pi 11 ) — 


■nr 3 N 2 

24 

nr 3 N 2 

6 


(8A 3 + 12mCA 2 + 6m 2 ( 2 A + m 3 ((4 + £ 2 )) , 
4 


A 3 + 3m 2 (A + m 3 ( 2 


2-C 


where F^ 1 ) and F( n ) are the free energies of phase I (A < A c ) and phase II (A > A c 


(3.20) 


respec¬ 


tively. Using (3.20) we find that the free energy has a discontinuity in its third derivative at 


the critical point 


d A (F (//) - F (/) ) 


A=A„ 


dl (F<"> - F«) | A . = 0, 


d l { F(II) - F {1) ) | A _A = 2ttt 3 N 


c 


C-2 


(3.21) 

(3.22) 


Thus at the critical point A c = (l — m the theory experience third-order phase transition. 
Third-order phase transitions are typical for matrix models in the planar limit and has been 
obtained for many different systems [34,261. As we have mentioned before there are examples 
of the phase transitions similar to the one we have described in this section. These are 


phase transitions in 3D Chern-Simons-Matter theory 24,25 and 4 D super-QCD with 2Nf 
massive hypermultiplets [2l], In both examples the phase transition of the theory in the 
decompactification limit at large -N is of third order. 



Figure 2: The free energy of AT = 1 SYM with fundamental hypermultiplets. The orange 
dots show the results for the numerical solution, while the dashed blue lines represent the 


analytical solution (3.20) 

. The parameters of the theory are taken to be r = 40, m = 6, ( = 3 , N = 100 corresponding 
to A r « 3.77. 


We have also checked results for the behavior of the free energy for different parameters 
numerically. In Fig. [2] we show the numerical and analytical results for the certain set of 
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the parameters. The orange dots show the results of the numerical simulation, while the 


dashed blue lines represent the analytical result (3.20). As we see the analytical results are 
consistent with the numerical ones. 


However from the Fig,2(b) we see that the third derivative, obtained numerically, is not 


exactly discontinuous but rather interpolates smoothly between two constant values, thus 
leading to the absence of the third-order phase transition. This happens because numerical 
calculations were performed for the large, but finite r. The finite volume of the system 
leads, as usually, to the absence of phase transitions, that are substituted by the crossover 
transition. However numerical study of the curve of d\F for different values of r have shown 
that it converges to the step function in the decompactihcation limit r —» oo, signaling about 
third order phase transition. 


3.4 Wilson Loops 

We also study the behavior of the circular Wilson loops near the critical point. In the planar 


limit the exponent term in the matrix model expectation value (2.18) does not affect the 


position of the saddle point, which leads to the following integral expression for the Wilson 
loop vev: 


(W) = / 


(3.23) 


Calculating this integral using the eigenvalue density p(0) from (3.13) and (3.16) we arrive 
to 


(W) 1 ' 1 ' 1 = cosh(27rra) , 


for the phase I and 


(W) (n] — ^1 — cosh(27rra) + cosh(27rrm) 


(3.24) 


(3.25) 


for the phase II. Using these expressions we find the discontinuity in the second derivative 
of the Wilson loop expectation value with respect to A at the critical point A c 


d\ - (H0<'>)L A ~ - - 4n 2 r 2 cosh(2;rrm). 

c 2 — ( 

which is consistent with the third-order phase transition. 


(3.26) 


4 J\f = 1 SYM with adjoint hypermultiplet 

In this section we consider phase structure of = 1 SYM with the massive adjoint hyper¬ 
multiplet of the mass m. As we will see this phase structure consists of infinite number of 
third-order phase transitions on the way from weak to strong coupling. 
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In the case of the adjoint hypermultiplet using (2.11) we obtain 


2 _ 2 
9 e ff ~ 9 ym > 


(4.1) 


i.e. in this case the coupling does not get shifts from the UV cut-off A and from now on we 
use the bare coupling constant Qym- 


In the planar limit matrix integral (]2.16|) is dominated by the saddle point 

dF(fa 


,, 7T ( ^ $ + S \ F v(4>i ~ h) + - h) 

dfa d(j>k \9 ym i 


= 0. 


(4.2) 


Using expressions (A.5) for the derivatives of Fy and Fh we obtain the following equations 
of motion 


167r 2 r 

9ym 


fa = ^ [f2 - (fa - 0j) 2 ) coth(7r(0j - fa))+ 


1 / 1 


- + (fa - fa + m) J tanh(7r(0j - fa + m))+ 

^ + (fa — fa — m) 2 ^ tanh(7r(0i - fa - m)) 


(4.3) 


2 V 4 


Restoring the radius dependence (2.19) and taking the decompactihcation limit, we obtain 
167T 2 , 


9ym 


E 


-(fa - 0j) 2 sign (fa - fa) + -(fa - fa + m) 2 sign(0j - fa + m)+ 


-(fa - 4>j - mf sign(0i - fa - m) 


(4.4) 


These equations can be directly observed from (2.20). Their continuous limit has the follow¬ 
ing form 


I67U 


0 = / di(jp(fa 


-(0 — 0) 2 sign(0 — 0) -|—(0 — 0 + m) 2 sign(0 — 0 + m) + 


— (0 — 0 — m) 2 sign(0 — 0 — m) 


(4.5) 


where we have introduced’t Hooft coupling constant t = g\ M N p| 

To solve (4.5) we proceed in the same way as while solving (3.8) for the Af = 1 SYM with 
the fundamental hypermultiplets. We take three derivatives of (4.5) , thus reducing integral 
equation to the algebraic one. Then we find general solutions to this algebraic equation for 
different phases of the theory, corresponding to different values of the parameters. Finally, 


5 Note that this’t Hooft coupling is different from A = g'y M N/r used in |6 
dimensionless constant, while t used in this paper has the dimension of length. 
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we fix all integration constants and free parameters by substituting general solution back 
into original equation. Three derivatives of (4.5) are given by 


167T 2 

t 


d'l/jpfy) [—2 o — ip\ + |0 — ip + m\ + |0 — ip — m |] , 


(4.6) 


0=J d'lJjpt'ijj) [—2 sign (0 — 0 ) + sign (0 — 0 + m) + sign (0 — 0 — m)\ , (4.7) 

0 = —4p(0) + 2p(0 + m) + 2p(0 — m) . (4.8) 

Before we start solving these equations let’s understand what do we expect to obtain. In 


analogy with 3D and 4 D results 26,19,21 


We start with the theory with 2 a < m, where a is the position of the support endpoint 


for the solution of (4.5). 


Then we increase the ’t Hooft coupling t which leads to the widening of the support 
(increase of a). At some value of t, where 2 a = m, we obtain emergence of two 
resonances, originating from the nearly massless hypermultiplets, i.e. from the terms 
|0 — 0| « m in the saddle point equations (4.5). Appearance of these resonances leads 


to the change of the form of the eigenvalue distribution, signaling about the emergence 
of a new phase. 

If we further increase coupling new secondary resonances appear every time when 
2 a = nm with n G Z. Each time this happens we expect to obtain phase transition in 
the theory. 

In the limit of very large ’t Hooft coupling t —> oo we expect to obtain an infinite 
number of resonances that will be averaged to the strong coupling solution obtained 
in ( 6 l 


Our final goal here is to obtain solutions of (4.5) for arbitrary number of resonances n. 


However in order to understand how to construct this general solution we start with finding 
explicit solutions for n = 0 and n = 1 and only then generalize solution to general n. 


4.1 Simple examples 


The phase with no resonances corresponds to the case of 2 a < m. To solve (4.5) with 


these values of the parameters we address (4.7) corresponding to the second derivative of the 


original equation. This equation then is equivalent to (3.12) and hence can be solved by 


p^(<t>) = -(S(4>-a)+5(4> + a)) 


(4.9) 


Here and further below we treat the h-functions at the endpoints of the distribution in the 
same manner as we were doing it in the case of the fundamental hypermultiplets, i.e. we 
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always assume the small shift in the argument of the 5-function that brings whole 5-function 
inside the eigenvalue support. 


Substituting ansatz (4.9) into the original equation (4.5), we obtain the following relation 


for the position of the support endpoint 


(o) 87f2 

a K ' = m - 

t 


(4,10) 


By construction this solution is valid when 2< m or, equivalently t < Thus the 

first critical point is 


t^ = 


167r 2 


m 


(4.11) 


Similarly to the case of the fundamental hypermultiplets to justify this solution we consider 


see appendix |D} 


solutions to the full equations of motion (4.3) in the limit of large but finite r, and then 
obtain (4.9) together with (4.10) in the limit of infinite r. For the details of this calculation 


We also found numerical solutions to (4.3). We show these solutions in Fig,3(a) with the 


orange dots. The dashed blue lines on the same plot represent positions of the 5-functions 
in our analytical solution (4.9). As we see from this picture analytical solutions reproduce 


numerical ones very well. All differences between these solutions appear due to the effects 
of finite r, which are discussed in appendix |D} 

Now we increase the ’t Hooft coupling and pass the critical point t^\ so that now 
m < 2a < 2m. Under this condition we expect two resonances to appear at 0 = ±(m — a). 
To solve equations of motion (4.5) for this case we divide the eigenvalue support [—a, a] into 
three regions 


Pa{(P), 

p (1) (0) = p B (0), 

Pc(0), 


m — a < 0 < a , 

—m + a < (f) < m — a, 
—a < (f) < —m + a , 


(4.12) 


and find general solution to (4.5) in every region. We also use 0 —> —0 symmetry of (4.5), 
that leads to 


pb(4>) = pb(-<!>) , pa(4>) = pc(-4>) 


(4.13) 


Region A (m — a < (f) < a). On this interval equation (4.8) reads 

-Ap A (<j>) + 2p c (<j>-m)=0, (4.14) 

where we used p(0 + m ) = 0 as 0 + m > a and also (0 — m) G C. Using the symmetry 
property (4.13) we obtain 


2 p A {4>) = pA{m - 0 ) • 


(4.15) 
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This equation can be satisfied only if Pa(4>) = 0. In analogy with the previously found 
solutions we assume that this is true up to isolated points on the boundary of the interval 
at (f> — a, rri — a. Then the natural ansatz is 

Pa{ 4>) — cofi(4> ~ a) + <q5(0 - m + a). (4.16) 


P(0) P<0) 




(a) t = 12 . (t < 


(b) t = 18. (t^ <t< 4 2) ) 


Figure 3: Free energy of the AT = 1 SYM with adjoint hypermultiplets. The orange dots show 
the results for the numerical solution, while the dashed blue lines on (a) and ( 6 ) represent 
positions of the 5-functions in the analytical solutions (4.9) and (4.27) respectively. The 
parameters of the theory are taken to be r = 15, m — 10, N — 100 corresponding to 
t { c ] « 15.79 and fl 2) « 23.68. 


Region B (—m + a < 0 < m — a). In this region (4.7) takes the form 

J dipp(ip)sign((j) — ip) = 

a—m 4> m—a a 

/ dijrpcW + / chpp B (i 0 - / di>p B (i’) - / <%uW0 = o 


(4.17) 


The first and the last terms in this expression cancel each other due to the symmetry prop¬ 
erties (4.13) of p(4>). The remaining part of the equation then reads 

(j) m—a 

[ dippsiip) = [ dipp B {i>), (4.18) 


suggesting, similarly to the case of n — 0 , that solution for the eigenvalue density P_b( 0 ) is 
sharply peaked around endpoints of the interval B at 0 = ±(m — a) 


Pb( 0 ) rs./ 5(0 — m + a) + 5(0 + m — a) 


(4.19) 
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Note that we have already took one of this 5-functions into account in Pa{4>) i n (4.16), while 
the second 5-function appears in pc(4>) as we will see further. 

Region C (—a < 0 < a — m). Finally, on the last interval C we do not need to solve 
any equations in order to find the eigenvalue density. Instead we utilize our previous result 
(4.16) for the density in the region A together with the symmetry properties (4.13) to obtain 

Pc(0) = Cq5(0 + a) + Ci<5(</> + m- a). (4.20) 


Notice here that, just as we expected, the second 5-function from p B (0) showed up here. 

Summing up all our results we find the following ansatz for the eigenvalue density in the 
phase with one pair of the resonances 

= c o 5(0 — a) + Ci5(</> — m + a) + Ci5(0 + m — a) + c o 5(0 + a). (4-21) 


This general solution has three free parameters cq, C\ and a. To fix all of them we use 
the normalization condition (3.14) together with the equation (4.6). The consistency of 
the obtained solution with the original equation (4.5) can then be checked by the direct 
substitution. From the normalization condition we obtain 


Co + Ci — - , 


(4.22) 


while from the equation (4.6) it follows 
167T 2 


t 


= 2 a (— 2 cq + ci) + 2 co(a + m) + 2 ci( 2 ?n — a) 


+ (—2d + c 0 ) (|0 - a + m\ + |0 + a 
This equation can be satisfied for any 0 only if 


m\) 


(4.23) 


Co — 2ci. 


(4.24) 


Combined with equation (4.22) this gives 


Co q j Ci 

3 6 


(4.25) 


Finally substituting these values back into (4.23) we obtain the following expression for the 
endpoint position 


= 2m — 


24vr 2 


(4.26) 


Corresponding eigenvalue density is given by 

P (1) (0) = ^(0 - a) + t)5(0 -m + a) + l5((j) + m-a) + ^5(0 + a). (4.27) 

3 6 6 3 
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By construction this solution is valid, when m < 2a < - 1 - ) < 2 m or equivalently < t < 

’ ^ J m m 

Thus the second critical point is given by 


(?> = 


24V 


m 


(4.28) 


As in the previous cases we justify this solution by finding the solution of (4.3) for large 


but finite r and taking r —> oo limit in the very end of the calculations. As the result we 


reproduce solution (4.27) with the right positions and the coefficients of the h-functions. For 


the details of this calculation see appendix [D] 


In Fig.3(b) we compare the analytical solution (4.27) with the corresponding numerical 
solution of (4.3) represented by the orange dots. The dashed blue lines show the positions 


of the h-functions in (4.27). As we see the analytical solution coincides with the numerical 


one up to the effects of finite r. 


4.2 General solution 

Now, using examples of the solution with n = 0,1 number of resonance pairs, we are ready 
to find a general solution with an arbitrary number n of the resonance pairs. This solution 
is valid when inequality mn <2 a < m{n + 1), where n € Z, is satisfied. In this case 
the cut can contain n pairs of the resonances. The mechanism of the emergence of these 
resonances is exactly the same as it was in the case of n = 1. Primary resonances appear at 
m) due to the discontinuities in the kernel of (4.5), caused by the peaks at the 


= ±(a 


endpoints of the eigenvalue distribution 0 = ±a. These resonances in turn create secondary 
resonances at 0 = ±(a — 2m) by the same mechanism and so on up to the last resonance 
pair at 0 = ±(a — nm) that fits inside the eigenvalue support. 


Thus generalizing solutions (4.9) and (4.21) to the arbitrary number of resonances we 


expect to obtain the eigenvalue density of the following form 


p( 0 ) = 5^ c k (£(0 — a + mk ) + h (0 + a — mk )) . 


(4.29) 


k =0 


Now we proceed in the same way as for the cases of n — 0,1. Namely we substitute ansatz 


(4.29) into equation (4.6). This substitution results in 


167T 2 


Cfc [—210 — a + mk\ + |0 — a + m(k + 1 )| + |0 — a + m(k — 1 )| + 


k =0 


n— 1 


n+1 


(4.30) 


(a ->■ -a, m ->• —m)\ = -2 ^ c k + ^ c k+ 1 + ^ c k -1 J (|0 — a + mk\ + |0 + a - mk\) , 


k =o 


k =-1 


k =1 


where in the last line we have shifted the summation index for the terms |0 ± a =F m{k + 1 )| 
and |0 ± a =F m{k — 1)|. Assuming that c_i = c n+ \ e 0 we arrive to the relatively compact 
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equation 


—— = y^C k {\(j)-a + mk\ + \(p + a-mk\) + 

t t^o ( 4 - 31 ) 

c n (10 - a + m{n + 1)| + |0 + a - m(n + 1)|) + c 0 (|0 — a — m\ + |0 + a + m) |) , 
where for simplicity we introduced 


Cfc — Cfc+1 2Cfc T Cfc_i . 


(4.32) 


Notice now that for the last two terms the following relations are satisfied everywhere on the 
support [—a, a] 


|0 — a + m[n + 1 )| + |0 + a — m(n + 1 )| = 2 (m(n + 1 ) — a), 
|0 + a + m\ + |0 — a — m\ = 2 [m + a). 


(4.33) 


However the terms in the summation in (4.31) depend on 0 and, in analogy with the pre¬ 


viously considered cases, we expect the coefficients of these terms to be zero. To prove this 
statement precisely we follow the same algorithm as in the case of n = 1. However, we 
should consider cases of even and odd number of the resonance pairs separately. Though 
they are similar and lead, as we will see soon, to the same result, there are some differences 
that should be briefly discussed. 


I-I- 


- i i V i i —h 


, 0 * 


N'T. 




ON 





yA o y 

I-f—1-- lv' I V f —1 1- 

P y yA on 


(a) Even n 


(b) Odd n 


Figure 4: Break of the eigenvalue support into intervals for the cases of even and odd number 
n of the resonance pairs. 


In the case of the even n the entire support [—a, a] should be separated into different 

There are 


regions with the boundaries corresponding to the resonances as shown in Fig 4(a) 


two types of intervals we should consider. Intervals of different types are shown with different 
colors in Fig. 4(a) On this picture and everywhere further for shortness we introduce the 


notation ak = a — mk. 
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Intervals < (ft < —a^ + j + i , j = 0,| — 1 (blue intervals in Fig,4). If (ft belongs to 
one of the regions we obtain the following relations 

, n , n 

20 , 7,-J <k<-+], 


n 


|(ft — a + mk | + \<ft + a — mk\ = 2 (a — mk ), k < — — j — 1, 


(4.34) 


n 


—2(a - mk), k > - + j + 1, 


Using (4.33) and (4.34) equation (4.31) can be rewritten as 

?+j 


I 67 H 

t 


=2 £ c» 

n „• 


2 3 ^ 


(4.35) 


+ 2 [ ^ ^ j Cj,(t7 — m/.'j + 2c,Jrn(n - li - a) + 2c (l i4; + m), 

^ — 0 fc= 7^ +J + 1 y 


Intervals —a^ + j < <ft < an_j , j — 1,.., | (green intervals in Fig,4). In these regions sim¬ 


ilarly to the previous case we obtain 


n , n 

20 , 2~j + l<k<-+j, 


n 


\<ft — a + mfc| + |0 + a — mk\ = 2(a — mk) , k < — — j 1 


(4.36) 


n 


—2(a - mk ), k > - + j + 1, 


and the corresponding equations of motion (4.31) turns into 


167t~ 

t 


-~2 ^ Cfc 0 

fc=f-i+i 

/f-t n 

+ 2 I J Cfc(a — mk) + 2 c n {m{n + 1) — a) + 2c 0 (a + m) 


(4.37) 


^k =0 fc=f+.?+l, 

Now the algorithm of defining all (n + 1) coefficients Cf- is the following. First condition 
comes from the normalization condition (3.14), that reads for the ansatz (4.29) as 


2 Y,Ck = 1 . 


(4.38) 


k=0 


We then need n more conditions to determine all the coefficients. These conditions can be 
obtained by considering the intervals shown in Fig, 4(a) We start from the region closest to 
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the origin, i.e. as < 0 < -a» + i, corresponding to j = 0 in (4.35). Then the first sum in 
this equation is given just by one term 


5+J 

2 y] Ca, 0 —> 2 Ch 
*=f-j 


(4.39) 


As before we assume that this ^-dependence should be canceled leading to Cs = 0. 

Then we move away from the origin to the next interval — a| + i < 0 < os_i. This interval 
corresponds to the equation (4.37) with j = 1, which turns the first sum into the sum of two 
terms 


\+3 


2 ^ C k (j) -)• 2 (C|+C| + i)0, 


(4.40) 


thus resulting in the condition + C« + i = 0. Repeating this procedure for all n intervals 
on the positive 0 half-axis we obtain n equations 


Y.Ck = « 


(4,41) 


k=1 


For the case of the odd n all considerations are analogous. As in the previous case we 
consider the division of the whole eigenvalue support into two types of regions, which is 


shown in Fig, 4(b) 


On the intervals a^±i_ i < 0 < — an±i +J , ,.7 = 1, (blue intervals in Fig 


(]4.3l|) can be written in the form 
^+3 

2 ^ C fc 0 


4) equation 


167T * 2 






n \ 


W E 

- £ 

C k (a — rri, 

y k =o 

k=^+jj 


intervals - 

CL n+1 | j < (f) 

^ CLn — 1 • 

2 J 


(4.42) 


71—1 


same equations can be rewritten as the following 
167T 2 


4) the 


^+j 


t 


2 Yi 

k-Hk-i 


2 J 


(4.43) 


2 I C k (a — mk) + 2c n (m{n + 1) — a) + 2c 0 (a + m) 


k =o 


fc=2i±3 +i 


24 














Now as in the case of even number of the resonance pairs we consider intervals and corre¬ 
sponding equations (4.42) and (4.43) one by one. 

To obtain n conditions for the coefficients Ck we start with the region —an +1 < 6 < a»-i. 

_ 2 2 

This interval corresponds to the equation (4.35) with j — 0 resulting in Cn+i = 0. Repeating 
procedure for every single interval on the positive half-axis we obtain the same equations 
(4.41) for the coefficients Ck- 


Now employing equations (4.38) and (4.41) let’s determine the coefficients cy. in our ansatz 


(4.29). To do this notice that the general solution for the recurrence relation 

c fe+ i - 2c fe + c fe _i = 0 , k — 1, ..,n. 


is given by 


Cfc = a + (3k 


(4.44) 


(4.45) 


The easiest way to see this, is to write the recurrence relation for large k in form of the 
differential equation c(k) = 0 with the general solution given by c(k) = a + /3k. 

To determine the coefficients a and (3 we use boundary condition c n+ i = 0. Substituting 
k = n + 1 into (4.45) we obtain 


a = — /3(n + 1) • 


Now substituting (4.45) into normalization condition (4.38) we find 


J> = ot(n + 1) + /3 


n{n + 1) 


k =o 


1 

2 


Finally combining this relation with (4.46) we obtain 

1 


a = 


P = 


(4.46) 


(4.47) 


(4.48) 


n + 2 ’ ' (n + l)(n + 2) ’ 

The last ingredient we need for the completeness of the solution is the position of the 
cut endpoint a. To determine it we should come back to equation (4.31). As we have shown 
Ck = 0 for k = 1, ..,n, so that corresponding terms do not contribute into equations. Using 
a and f3 determined in (4.48) we obtain expression for Co 

1 


C 0 = Ci — 2 c 0 = (3 — a = 


n + 1 


Substituting it into (4.31) we find the position of the support endpoint 

47T 2 

a ^ = (n + l)77i- —(77 + 1)(?7 + 2). 

Obtained solution is valid, when mn < 2 < m(n + 1) or equivalently when 

8tt 2 8tt 2 

- (77 + 1) < t < - (n + 2), 

777 777 


(4.49) 


(4.50) 


(4.51) 
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Figure 5: The free energy of J\T = 1 SYM with adjoint hypermultiplets. The orange dots show 
the results for the numerical solution, while the dashed blue lines on (a) and ( b ) represent 


positions of the 5-functions in the analytical solutions (4.53). The parameters of the theory 
are taken to be r = 25, m = 5, N = 200. 


so that the n th critical point is located at 


4 n) = —(n+l) 
m 


(4.52) 


Let us summarize the findings of this section. We have obtained the general solution 

If the ’t Hooft couping is in the interval 


to the matrix model saddle-point equation (4.5) 

^2 , . 


+ 1) < t < + 2) the solution contains n pairs of resonances. The eigenvalue 

density for this solution is given by 


_^ ^ 

p{4>) = Ck W - a + rnk ) + 5((/> + a - mk )) , c k = -————— , 

(n + l)(n + 2) 


(4.53) 


with the position of the eigenvalue support endpoint a given by (4.50). 


If we take n — 0,1 in this general solution we will reproduce corresponding formulas (4.9), 


and (4.27) obtained previously. We have also checked general solution (4.53) numerically. 
In Fig|5|we show with the orange dots numerical solution to (4.3) for the cases of n = 3 and 
n = 5. On the same plots we show the positions of the 5-functions in our analytical solution 


(4.53) with the vertical dashed lines. As we see numerical and analytical results are in good 
agreement. The positions of the 5-functions reproduces positions of the peaks very well. The 
height of the resonances also decreases when moving away from the endpoint, as predicted 


in (4.53). 
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4.3 Free energy and the order of the phase transition 


Now we are ready to address the question of the order of the phase transitions taking place at 
the critical points For this we evaluate the free energy corresponding to the eigenvalue 


density (4.53). We read the free energy in the decompactihcation limit directly from (2.20) 


choosing appropriate matter content. In the continuous limit it can be written as follows 


F = 


1 


47T 


2nr 3 N 2 



F = — + — / / d0d0p(0)p(0) [10 - 01 


12 

~10 - 0 + m\ 3 — ^|0 — 0 — rn\ 3 


(4.54) 


where we have also introduced F, corresponding to the free energy with excluded dependence 
on the radius r and rank N of the gauge group. 


Integrals in the expression (4.54) can be evaluated after substituting our solution (4.53). 


The first integral is relatively easy to evaluate 

/<->w=^ i> 

J u —n 


— mk ) 2 = 


k =o 


2tt 


(4.55) 


—— (1 + n)(2 + n) [967r 4 (l + n)(2 + n) — 167r 2 fm(3 + 2 n) + 3m 2 f 2 ] , 
3 t A 


where in the last step we hav e used explicit expressions for the coefficients ry c (4.53) and the 
endpoint position (4.50). 


The second integral coming from the one-loop determinants is more complicated. Let’s 
consider it in more details here 


4'° = ~24 / P(0)P(0) 


-2|0 - 'ifj | 3 + 4> - 0 + m\ 3 + ^|0 - 0 - m\ 3 


12 


EE CkCi [—2m 3 \k — l \ 3 — 212a — m(k + Z)| 3 + m 3 \k — l + 1| 


(4.56) 


fc=0 «=o 


+?n 3 \k — l — 1| + \2a — m(k + l + 1)| 3 + |2a — m(k + l — 1)| 3 ] . 


In the same manner as it was done when solving the saddle-point equation in (4.30), we shift 
the summation index l in the last four terms to obtain 


4 n) = - ^2 CkCl i m3 \ k ~ l \ 3 + l 2a _ m ( k + 0I 3 ] ~Y2^2 Ck t C ° + ^ (a 

z k= 0 1=0 z fc=0 (4-b<) 

+ (2a — m(k — l)) 3 ) + c n (m 3 (n + 1 — k) 3 + (m(A: + n + 1) — 2a) 3 )] , 


where we have also used inequalities mn < a < m(n + 1). This expression is much easier 
to work with, because, as we know, Ck = 0 for k = 1, ..,n. Then substituting values of the 
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coefficients c 0 = c n = ( n+ y ^ n+2 ) , = — 7^1 and the endpoint position a^ from (4.50) 

we finally obtain 


An) __ 

2 12f 3 


[512(1 + n) 2 (2 + n) 2 7r 6 — 64mt7r 4 (l + n)(2 + n)(3 + 2 n) + m 3 f 3 ( 3 + 2n)](4.58) 


The free energy of the system is then given by the sum of (4.55) and (4.58) 


F^ [2567t 6 (1 + n) 2 (2 + n) 2 — 647r 4 mf(l + n)(2 + n)(3 + 2n) + 


247r 2 f 2 m 2 (l + n){2 + n) — m 3 f 3 ( 3 + 2n)] 


(4.59) 


To hnd the order of the phase transition we calculate derivatives of the free energy (4.59) 
with respect to the’t Hooft coupling at the critical points (4.52). After short computation 
we observe 


9i 3 (JTC+I) _ F <»>) | ( _ ( ,„ +1) 


a 2( F („+l)_ F W)|^ n+i ) =0j 

m 6 

_ 512tt 6 (2 + n) 3 ' 


(4.60) 

(4.61) 


Hence, when ’t Hooft coupling t is increased theory undergoes through an infinite chain of 
phase transitions of the third order. Moreover, with the increase of the coupling transition 
becomes weaker transforming into crossover transitions or infinite’t Hooft coupling t. Similar 
behavior was obtained in the case of the mass-deformed ABJM theory 26,27 . 



Figure 6: The free energy of J\T = 1 SYM with adjoint hypers. On the pictures above 
,_F , where F is the free energy defined in the usual way. The orange dots show 


1 


p — ___ 

2nr 3 N 2 J 

the results for the numerical solution, while the dashed blue lines represent the analytical 
solution (4.59) 

. The parameters of the theory are taken to be r = 25, m — 10, N — 100. 


























We have also evaluated the free energy for the full matrix model (2.16) with massive 
adjoint hypermltiplet numerically. The results of these numerical simulations are shown in 
Figj6]with the orange dots. The dashed lines on the same plots represent our analytical result 
( 4.59[ ) for the free energy. As we see numerical and analytical results are in good agreement. 
The only difference can be obtained in the graph of the third derivative of the free energy. 
As we see the free energy, obtained numerically, does not have discontinuities at the critical 
points and hence leads to the absence of the phase transitions. This result is expected 
because the numerical solution is found for the finite r and, as known, finite-dimensional 
systems do not experience phase transitions. 


4.4 Wilson Loops 

We also study the behavior of the circular Wilson loop near the phase transition points. In 
the planar limit contribution of the e 2n '^ prefactor in (2.18) doesn’t effect the position of the 


saddle point. Thus we express Wilson loop expectation value in the integral form 

(W) = f 


(4.62) 


Substituting solution (4.53) into this expression we easily obtain 


(IT ) (n> = 2 ^ Ck cosh ( 27r 


k =0 


47r 2 

(n + 1 — k)m -— (n + l)( n + 2) 


(4.63) 


where we have used expression (4.50) for the position of the endpoint. Calculating derivatives 
of (ITwith respect to the coupling constant around the critical point we obtain 

M<w0 (n+1) - <wo (n) ) Lh" +1 ) = ° 


d 2 t ((W) (n+1) - (W) (n) ) \ t=t (n+D = 


m 4 sinh ( irm(2 + n)) 


(4.64) 

167t 2 (2 + n) 2 sinh(7rm) 

Hence, the first discontinuity appears in the second derivative of the Wilson loop. This result 
is consistent with the third-order phase transitions experienced by the system at the critical 
points ■ 


4.5 The strong coupling limit 


It is interesting to study the behavior of the solutions (4.53) for a very strong coupling 
constant, such that 2 a m. In terms of our solution, this means that the eigenvalue 
density contains large number of resonance pairs 1. 

Previously in |7| we have obtained that for the strong coupling eigenvalue density at the 
saddle point is 


p{4>) = 


32vr 2 


(9 + 4?h 2 ) A ’ 
= 0, 


< a , 


> a 


(4.65) 
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where the endpoint position is given by a = (9 + 4m 2 ) A/647T 2 and A is the coupling constant 
defined as A = g\ M N/r = t/r. Restoring r-dependence (2.19) and taking the decompactifi- 
cation limit r —)■ oo we obtain 


P(0) = 


8vr 2 


m 2 t ’ 
= 0, 


< a , 
> a , 


(4.66) 


where the cut endpoint position is given by 


a = 


m 2 t 

167T 2 


(4.67) 


Now let’s compare the expressions above with the corresponding limits of our solution 


(4.53) and ( 4.50 ). In order to do this recall that the solution with n pairs of resonances is 
valid when ^-(n + 1) < t < —(n + 2). Thus for very large n we can approximate the 
number of the resonance pairs by n ps |^|. For the position of the endpoint using (4.50) we 
can easily obtain 


An 2 2 

a ~ mn - n 

t 


m 2 t 
167T 2 ’ 


(4.68) 


which exactly reproduces the expected result (4.67). 


In order to reproduce the eigenvalue density (4.66) we consider the case of large n in 


(4.53). As n 1 the cut contains a large number of the ^-functions which is then regarded 
as a continuous distribution. In order to understand which distribution is this, let’s consider 
the following integral 


d0p (n) (0) = Ck 


n 


dk —-— =-= 


n* 


x 

mn 


8n 2 x 

m 2 t 


(4.69) 


2 m 


where we have used (4.53) for p^ n \(j)) and C&, and approximated the sum over k by the 


integral assuming that the sum goes to large enough values of k. As we obtained linear 

X 

dependence of the integral f p(0), we conclude that the eigenvalue density is given by the 

o 

constant on the interval [—a, a] 


PW = 


87r 2 
m 2 t ’ 


= 0, 


< a, 
> a, 


(4.70) 


which exactly reproduces solution (4.67). 


To conclude this section, we have shown that for the large number of resonances we 
can consider our solution as the averaging over these resonances. This averaging results in 
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Figure 7: Results of the numerical solution (orange dots) to the equations of motion (4.3) 
with the parameters m = 5, r = 25, N = 200 and very large’t Hooft coupling. As we see 
for the large t solution approaches distribution (4.66), which is shown by dashed line 


the constant eigenvalue distribution obtained in |7|. Similar behavior was obtained for the 
solutions containing resonances in the mass-deformed ABJM theory 26,27 and 4 D M = 2* 


theory 19 


The same conclusion about behavior of (4.53) at very strong coupling can be drawn 
numerically. On the plots in Fig j7] the orange dots represent numerical solutions to (4.3). On 
the same plots we show the solution ( ]4.66 ) with the dashed blue line. As we see distribution 
for t = 300 still looks like a bunch of peaks, while for t = 750 it already starts approaching 
(4.66) and for t = 1250 completely coincides with it up to small regions of the size ~ ^ around 
the endpoints of the distribution. Hence, numerical simulations support the conclusions of 
this section. 
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5 Squashing spheres 


Another issue to discuss is the effect of the geometry on the phase transitions described in 
this paper. Localization methods can be generalized to theories on the space with geometry 
more complicated than the spheres. In particular 5 D SYM was considered on squashed 

Here we focus on the squashed spheres, which can be embedded into C 3 as follows 


spheres 35,36 , Y p,q spaces 19,10 and even more general toric Sasaki-Einstein manifolds 11 


2 I 12 , 2 I 12 I 2 \ \2 2 

UJl | Z\ | + cu 2 \z 2 \ Y w 3 | Z 3 \ — V . 


(5.1) 


Then the corresponding matrix model obtained after localizing 5 D SYM on this space is 

(5.2) 


f -^flipTr( 0 2 )-YpTr( 0 3 ) / \ , 

Z = / d(f)e 9 ym detndf S 3 (i<j) uj)) det R [S 3 


3 

+ 2 


u 


where u = (cni, uj 2 , oj 3 ) and p is the volume factor of the squashed sphere, which is given 
by the ratio of the volumes of squashed and usual spheres 


V°1 3%. 1 

Vol S 5 i>Ji oj 2 UJ 3 


(5.3) 


Finally, S 3 denotes triple-sine function that can be defined through the following infinite 
product 

OO 

S3 (x| uj) = (ji\ + n 2 UJ 2 + n 3 a) 3 ) (('Ui + 1) uq + (n 2 + 1) uj 2 + ('n 3 + 1) cu 3 — x) .(5.4) 

111,2.3=0 


The matrix model (5.2) was partially studied for the case of 5 D SCFT, i.e. SYM theory with 
the USp(N) gauge group and infinite coupling constant (Jym 15,16 . One of the reasons to 


study supersymmetric theories on the squashed spheres is that their partition functions are 
related to the supersymmetric Renyi entropy. The latter can be evaluated from the partition 
function of theory on the n-covering of sphere, which is equivalent to the squashed sphere 
with the squashing parameters (4,1,1) 


17 


In general the matrix model (5.2) is very complicated and can not be solved directly. 
However, we are only interested in the behavior of the matrix model solutions in the decom- 
pactihcation limit r —>■ 00 which implies that the arguments of the triple-sine function are 
large, 1. Thus, we can use asymptotic expression for these functions 


16 


36 


log S 3 (ix\ o;)|, 


7T 


-x 3 - 


X 7T CJtot. 2 
-X 


sign(x) , 

_ 6cni u 2 UJ3 4o;i 0J2 uj 3 

77 + u lA+ <^1 <^3 + u 2 ^ 3 ) . 1 + U 2 + w i) 


I2uj\ uj 2 0 J 3 


-x + 


24uji uj 2 uj 3 


+ o(- 

x 


(5.5) 
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where u ta t. = <^i + ^2 + <^ 3 - In the case of the matrix model (5.2), x = 0 = r<p and we keep 
only cubic terms so that the partition function can be rewritten in the following form 


Z = 




OJI CO 2 C0 3 

2n r 3 


Cartan 

4 TT 2 


(5.6) 


F (</>) = ^2 tr (0 2 ) + —tr Ad |0| 3 - — tr «/ (10 + m \ 3 + 10 ~ m ! 3 ) 


^e// 


24 


From this expression we can see that the squashing does not affect our conclusions about 
phase transitions neither for fundamental nor for adjoint hypermultiplets. The only contri¬ 
bution of the squashing comes from the volume prefactor of the free energy. Similar results 
were obtained in the strong coupling limit of 5 D Af = 1 SYM with adjoint hypermultiplet 
on Y p,q spaces [9| and toric Sasaki-Einstein manifolds [II]. In particular, it was shown that 
for both cases the only difference in the prepotential of the theory between these spaces and 
five-sphere is in the overall volume factor. 

Note that effects of the squashing on phase transitions were also studied for the case of 
AD AT = 2* SYM in 23 . However results obtained in that paper differs from ours. The 


main conclusion is that the squashing of the sphere shifts all critical points of the phase 
transitions chain to the smaller coupling constant. While in our case as we have seen above 
squashing does not affect positions of the critical points. 


6 Discussion 

In this paper we studied the decompactikcation limit of the matrix model obtained from 5 D 
Af = 1 SYM theory with different hypermultiplet content on S 5 . First we considered Nf 
massive fundamental hypermultiplets coupled to Af = 1 SYM theory. Solving the planar 
limit of the corresponding matrix model we found that the eigenvalue density p(0) at the 
saddle point consists of an isolated peaks. These peaks are located at the boundaries of the 
distribution at 0 = ±a and at the points 0 = ±m. Hence, we concluded that there are two 
phases of the theory. One phase corresponds to the case when the support includes 0 = ±m 
and another one when it does not. Transition between these two phases takes place at the 
critical point 


t c = 


8 v^/ wy 1 
m V 2 NJ ’ 


( 6 . 1 ) 


where t = gy M N is the’t Hooft coupling used in this paper. Evaluating the free energy we 
find that the transition between these two phases is a third-order phase transition. 

Second theory we consider is Af = 1 SYM coupled to the massive adjoint hypermultiplet. 
In this case we also find that solution for the eigenvalue density consists of isolated peaks. 
These peaks are located at 0 = ±(a — mn ), where a is the endpoint of the distribution 
support, m is the mass of the hypermultiplet and n is an integer number. As the coupling 
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is increased the length of the support grows and it can include more and more peaks. As 
the result, theory goes through an infinite number of phases which differ by the number of 
peaks in the solution. Transition between phases with n and n + 1 pairs of peaks take place 
at 

8-7T 2 

j(») = _(„ + 1), (6.2) 

m 


Free energy calculation showed that these transitions are third-order phase transitions. In 
the limit of very strong coupling the peaks, forming the eigenvalue density, smooth out and 


solution becomes equal to the strong coupling solution derived in 18 


It has been also shown that the effect of the five-spheres squashing comes in the form of 
an overall prefactor of the free energy in the decompactihcation limit of the matrix model, 
while saddle point remains the same. Hence, we concluded that the squashing does not affect 
neither the positions of the critical points nor the order of the phase transitions. 

Finally, we note that the phase structure described in this paper seems to be general 
for the supersymmetric theories with massive hypermultiplets, as they arise in D = 3,4, 5 
dimensions. Hence, in the future it would be interesting to understand the nature of these 
phase transitions better. 
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A Properties of the functions f(x) and l(x 

In this appendix we consider properties of the functions /(x) and /(x) that appear in the 
expressions for the regularized one-loop contributions (2.12) and (2.13). Function /(x),first 
introduced in (3], can be written in th following form 


f(x) = 


mx 


x 


log (1 


r 2 ™) + —Li 2 (e- 2 ™) + -^Li 3 (e~ 2nix ) C(3) 
' 7T v ' An z v ' 


2vr 2 


while /(x), first introduced in |3l], is defined by 


l(x) = -xlog(l-e 2mx ) + - 


7 tx z H—Li 2 (e 

7T 


2ivix 


; 1 12 


Using these expressions we can derive the following properties of l(x ) and /(x) 
1 . /(x) is an odd function and /(x) is an even function 

l(x) = -Z(-x), /(x) = /(-x) 


(A.l) 


(A.2) 
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2. The derivatives of the functions are given by 


df{x) 


= nx 2 cot (Ax); 


dl(x) 


= —7 TX COtfAx) ; 


dx dx 

3. The asymptotic behavior of the functions is given by 

( j \ nr nr 

— V ix = -Ixl 3 H—Ixl; lim Imf (\ ± ix) = ±^x 2 ; 

2 J 3 4 x^roc V2 ’ 2 

lim Re l ( - + ix ) = — — |x|; lim Im l (4 ± i x \ = =p^x 2 ; 

|a:|—>oo V 2 J 2 x^roo V2 2 2 


7T , 


lim Re/ (ix) = — — |x| 3 ; 
|a"|—>oo 3 


Im/ (ix) = 0 . 


(A.3) 


(A.4) 


Using these properties we derive similar properties for F v (x) and F H (x) functions in (2.14) 


and (2.15): 


1. Both Fy(x) and Fh(x) are symmetric functions 

F v (-x) = F v (x), F h (x) = F h (x). 


2. The derivatives of the functions are given by 
dF v (x) 


dx 

dF H (x) 

dx 


= —— (2 — x 2 ) coth(7rx), 

= —— + (m + x) 2 ^ tanh(7r(x + m)) 


(A.5) 


3. Finally the asymptotic behavior of these functions is given by 


lim Fy(x) 

|:r|—>oo 

lim Fh(x) 

|ir |—>oo 


7T . . 

6 N 
7T , 


7T X 


7T , 


-x+m -x+m 

6 8 


(A.6) 


B Details of numerical analysis 

To obtain numerical results presented in this paper we have used the method introduced 


m 


37 in order to solve the ABJM matrix model. It was also adopted for solving 5 D SYM 


matrix model j?j and 5 D Chern-Simons matrix model 18 


Instead of solving the the saddle point equation = 0, where F is the prepotential, 
we introduce dependence of the eigenvalues on the effective ’’time” variable t. Then we 
substitute our original equation with the effective ’’heat” equation 


d(j)i(t) 

dt 


dF 

~dt 


(B.l) 
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At large time t solution of this differential equation approaches the solution of the original 
saddle point equation, provided the choice of r’s sign is made properly. 


However, solving (3.5) which corresponds to the case of the fundamental hypermultiplet, 
we faced some problems. Onr numerical method appeared to be unstable and very sensitive 
towards initial conditions due to the repulsive central potential ~ —A0 2 . 


To avoid this problem we used the fact that the renormalization (2.11) of the coupling 


constant ()ym for the case of J\f = 1 SYM can be reproduced by considering the case of one 
adjoint hypermultiplet with very large mass. The easiest way to realize this is to consider 


(4.3). If the mass of the hypermultiplet is large we can rewrite this equation in the following 


form 


167r 2 r 

9 ym 


2 mN ) 


(0i - 0,-) 2 j coth(7r(<^i - <pj)), 


(B.2) 


where we used tanh (7r(0* — <f>j ±rh)j = sign (±mir) = ±1 and also Yh&j = 0 due to the 

' _ ' j 

tracelessness condition. From (B.2) we conclude that the coefficient in front of 0, on the 


l.h.s. plays the role of effective coupling. To put (B.2) in the form similar to (3.5) A should 


be expressed through the bare coupling (]ym and the mass of the hypermultiplet m in the 
following way 

87T 2 

A = 77i - . (B.3) 

where t = gy M N is the’t Hooft coupling. So in order to avoid the problems with the repulsive 
central potential we introduce adjoint hypermultiplet with large enough mass m. In this case 


the central potential itself is attractive as gy M > 0. However the effective coupling (B.3) can 


be made negative as we need it in equation (3.5). All numerical simulations of the theory 


with the fundamental hypermultiplets presented in this paper were performed in this way. 

C Solution for the finite r: fundamental matter 


In this appendix we consider the solution of (3.5) for the large enough but finite r. In section 


[3] we have neglected all terms subleading in 1/r in this equation in order to consider properly 
the decompactification limit r —y 00 . As the result we have obtained solutions (3.13) and 


( |3. 16 ) containing 5-functions. These 5-functions arise in the solutions because of the kernel 
in (3.8) that looks like — (0 —0) 2 sign (0 — 0) leading to the attractive potential between 


the eigenvalues. If not the repulsive central potential, the eigenvalues would just all condense 
at the origin 0 = 0. 

Here we will be more accurate with the kernel of integral equation. We now do not 


neglect the subleading term in the kernel of (3.5), which leads to 


—2Ar 2 0 = -£r 2 (0 + m ) 2 sign(0 + m ) + -Cr 2 (0 — m) 2 sign(0 — m) + 


d0p(0) (2 — r 2 (0 — 0) 2 ) sign (0 — 0), 


(C.l) 
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where p(0) is the eigenvalue density defined by (3.7). Due to the symmetry of (C.l) we expect 
the eigenvalue density to be symmetric. Note that terms subleading in 1/r are repulsive, 
hence they wash out 5-functions into peaks of the finite width ~ 1/r. At the same time 
we have not included terms subleading in 1/r into the central potential, because they add 
small central repulsion, thus never playing an important role. Assuming cotli (7rr(0j — 0j)) « 
sign(0j — (j>j) we also neglect terms exponentially suppressed in large r. However, note that 
this approximation works pretty well even in the region, where \(f) t — 00 ~ 1/r due to the 
factor of 7 i in the argument of coth. Now we consider phases I and II separately. 


C.l Below the transition point 


We start with the phase below the transition point which corresponds to a < m , where 
a is the position of the eigenvalue support endpoint. In this case equation (C.l) and its 
derivatives are given by 


(2A + £ nn) r 2 <j) + j (2 — r 2 (0 — -0) 2 ) sign (0 — -0) 

= o, 

(C.2) 


= o, 

(C.3) 

4p'(0) — 2r 2 J d0p(0)sign (0 — 0) 

= 0, 

(C.4) 

4/(0) - 4r 2 p(0) 

= 0, 

(C.5) 

General symmetric solution to (C.5) is 



p(0) = C cosh(r0). 


(C.6) 


Now we are left to determine constant C and the position of the support endpoint a. To de¬ 
termine these constants we substitute the general solution (C.6) into normalization condition 
(p4| 


a 

j d(pp{(j )) 

—a 


-C sinh(ra) = 1 
r 


(C.7) 


and into equation (C.3) 


(2A + (m ) r 2 — 2 ar 2 + 4 C cosh(ra) = 0 . 


This leads to 


C = 


2 sinh(ra) ’ 

(2A + (m ) r 2 — 2ar 2 + 2r coth(ra) = 0 . 


(C.8) 


(C.9) 
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The consistency of these conditions with equations (C.2) and (C.4) can be checked by the 
direct substitution. In general (C.9) can be solved only numerically for different values of 
the parameters r,m,( and A. However, for the large radius, when ra 1, we can simplify 
the last equation 


a = A + -(m, 


(C.10) 


which reproduces the support endpoint (3.15) we found in section [3j Finally, the eigenvalue 
density below the phase transition point A c |^| is 


P {1 \4>) = 


2 sinh(ra) 
= 0, 


cosh(r</>), 


< a 
> a. 


(C.ll) 


For the large r this distribution is the function with two sharp peaks of width ~ 1/r at 
(j) ± a. In the decompactification limit r —s- oo this density is equivalent to (3.13), because 
the width of the peaks goes to zero, while the integral of p(0) stays 1. 

with 


We have also found the numerical solutions to (3.5), which are shown in Fig,8(a 


the orange dots. Dashed blue lines shows the solution (C.ll) with the coefficient C and the 
endpoint a found numerically from the equations (C.9). As we see, solutions found above 
are consistent with the numerical results. 


p(4>) 



p(4>) 

1.2 




a=3.44 — aw—7.5 -6 —m 



(b) A = 5.5 (A > A c ) 


Figure 8: The eigenvalue density p((p) calculated with the following value of parameters: 
r — 5, m = 5, C — |, N — 100 corresponding to A c « 3.76. Orange dots show the results for 
the numerical solution, dashed blue lines show the solutions (C.ll) and (C.25) on pictures 
(a) and (6) correspondingly. 


6 Notice that this critical coupling is in general different from (3.181 and can be found numerically using 


(C.9). However, in the decompactification limit r —> oo the critical point A c approaches (3.18). 
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C.2 Above the transition point 


Now we find the solution to (C.l) for the case of the system in the phase II. This phase 
corresponds to a > m, and to solve equations of motion we should divide the whole eigenvalue 
support [—a, a] into three parts A, B, C 


Pa(0), 

m < 0 < a ; 


p (//) (0) = p B (0), 

— 77i < 0 < m ; 

(C.12) 

Pc(0), 

—a < 0 < — m ; 


with the symmetry constraints 

PbO) = Pb{~4>) , 

Pa(0) = Pc{~4>) • 

(C.13) 


Now we write down equation of motion (C.l) in regions A and B and solve them in this 
regions separately. 

Region A (m < 0 < a). In this interval (C.l) and it’s first, second and third derivatives 
are 


2Ar 2 0 + ^<yr 2 (0 2 + m 2 ) + J d'ippfy) (2 — r 2 (0 — 0) 2 ) sign (0 — -0) 

2Ar 2 + ( r 2 0 + 4p^(0) — 2r 2 j dipp{tji) |0 — 0| 

(r 2 + 4p' A (0) — 2r 2 J dxppfy )sign (0 — 0) 

VA(0) - 4r 2 p^(0) 


0, (C.14) 

0, (C.15) 

0, (C.16) 

0. (C.17) 


As we do not have any symmetry restrictions on the general solution to 
by 


(C.17) 


is given 


Pa{4>) = Ci cosh 



a + 7n\ 



+ C 3 sinli 


^r(j) — r 


a + m\ 

) 


(C.18) 


Note that we choose solution to be centered around 0 = (a + m)/ 2. This choice is made for 
convenience. The shifts in cosh and sinh can be adsorbed into constants and C 3 

Region B (— m < 0 < m). On this interval 0 < m, so that everything is the same as in 
the case of the phase I. Equation of motion and its three derivatives then can be written in 
the form of (C.2)-(C.5). Similarly to the phase I we obtain 


Pb{4>) = Ci cosh(r0), 


(C.19) 


Region C (—a < 0 < —m). In this region using symmetry properties (C.13) we find 

a + m\ _ . , / a + m 


Pc(4>) = C 2 cosh r0 + r 


C 3 sinh 


+ r- 


(C.20) 
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To complete our solution, we need to find all the parameters C\, C 2 , C 3 , a in (C.18),(C.19) 
and (C.20). To determine these parameters we substitute solutions (C. 18)-(C.20) into equa¬ 
tions (C.2)-(C.4), (C. 14)- (C.16) and normalization condition (3.14). This leads to 

(C.21) 
(2-C )r, (C.22) 


Ci sinh (rm) + 2C 2 sinh ( r 


— yyi 

46*2 sinh [ r —-— ) + 4C 3 cosh I r 


a — m 
2 

a — m 


r 

2 ’ 


4 (C 2 — arC 3 ) cosh ( r 


a — in 


+ 4 (C 3 — arC 2 ) sinh ( r 


Q - YYl 

Ci cosh(mr) — C 2 cosh ( r - ) + C 3 sinh ( r 


2 

a — m 
2 

a — m 


- 2 Ar 2 , (C.23) 

0, (C.24) 


where (C. 21 ) is obtained from the normalization condition, while (C.22),(C.23) and (C.24) 
are obtained after substituting the solution for the eigenvalue density into (C.16), (C.15) 
and (C.3) correspondingly. It can be checked by the direct calculation that remained equa¬ 
tions (C.2),(C.14),(C.4) will be satisfied automatically if the coefficients C \, C 2 , C 3 , a obey 
equations we have found above. 

To summarize we have found that above the phase transition point, when m < a, the 
eigenvalue density satisfying equation of motion (C.l) can be written in the form 

a + m' 


C 2 cosh I rep — r - ) + C 3 sinh [ reft — r- 


pCT( 0 ) = Ci cosh(r 0 ) 


C 2 cosh ref) + r 


a + m 


— C 3 sinh ref) + r 


a + m 


77i < ep < a ; 

- m < ef> < 77i] (C.25) 

—a < ef> < —m ; 


with coefficients Ci, C 2 , C 3 and the endpoint a of the cut satisfying equations (C .21 )-(C.24). 

Equations (C.21)-(C.24) can be solved exactly in the limit of large radius r. In the case, 
when rm 1 and r^C^ > 1 we approximate these equations with simpler ones 


Ci exp (mr) + 2C 2 exp 
(C 2 + C 3 ) exp 
(C 2 + C 3 ) ar exp 
Ci exp (mr) + (C 3 - C 2 ) exp 



r, 



2A r 2 , 

0 . 


(C.26) 
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This system can be solved leading to 

a = A , c i = ^C^exp(-mr), 


= r|'i-ic)exp(r=^ 


Note that to find these solution we assumed 


C* = r - 


3 \ ( ttt — a 

-C exp r 


. (C.27) 


r a^n ^ x 


which means that solutions only 
« m. However this approximation 


makes sense far enough from the critical point where a 
will work very well in the decompactification limit even for the points close to the critical 
ones. 


As we see from (C.27) the position of the cut endpoint is consistent with the one obtained 


in the sectio n [3] I f we now come back to the eigenvalue density (3.13) and take into account 
coefficients (C.27) we can see that in the decompactification limit r —y oo solutions are 


nonzero only at the points 0 = ±a and 0 = ±m, where they take infinite value. At the same 


time from the normalization condition (3.14) we know that p(0) integrates to one. Then the 
natural function to describe this limit is (3.16) and the coefficients in front of the (5-functions 
should be 

p(0) = (C *2 + C 3 ) <5(0 — a) + (Ci + C 2 — C 3 ) 5(0 — tit) + (a —> — a, ttt —> —ttt) (C.28) 


which leads exactly to (3.16). 


We also can find solutions to (C.21)-(C.24) numerically for the particular combinations of 


the matrix model parameters. I 11 Fig, 8 (b) the orange dots shows numerical solution to (3.5), 


while the dashed blue lines shows the solution (C.25) with all coefficients found numerically 


from (C.21)-(C.24). As we can see numerical solution perfectly coincides with the solution 


found in this section. 


D Solution for the finite r: adjoint matter 


In this appendix we consider solution to the saddle-point equation of the matrix model (2.16) 


for the large but finite r. We have solved these equations in the decompactification limit 
(see section [4]) neglecting all terms subleading in 1 /r. As the result we have obtained the 
general solution ( |4.53 ) containing the 5-functions. In order to show more precise how these 
5-functions arise we should take into account first subleading terms in the kernels of (3.5), 
which leads us to 


1670 


r 2 0 = / d 0 p( 0 ) [(2 — r 2 (0 — 0 ) 2 ) sign (0 — - 0 ) 


- ( —b r 2 (0 — 0 — ttt ) 2 ) sign (0 — -0 — ttt) + (ttt —> —ttt) 


(D.l) 


where we restore the radius dependence (2.19) and assume that coth( 7 rr( 0 —0)) ~ sign(0—0) 
as well as tanh( 7 rr (0 — 0 ± ttt)) « sign(0 — 0 ± ttt) when we consider the limit r>l. The 


same equation can be obtained by a minimization of the free energy ( 2 . 20 ). 


41 




































Note that subleading terms correspond to the repulsive force between the eigenvalues at 
~ 1/r distances. Hence these terms wash out 5-functions into the peaks of 1/r width. As 
we will see further, it is possible to obtain the analytic form of these peaks and reproduce 


(4.53) in the decompactihcation limit r —> oo. 


To obtain solutions of (D.l) we first reduce this integral equation to the differential one 
taking three derivatives 


167T 2 4 1 , . 1 . . I 

— = ^PW) + ^PW + m ) + -fr2PW^ m ) +J 

+ |0 — 0 + m| + |0 — 0 — m\] , 

0 = ~ 2 m + + m) + p '(0 - m) + 

/j^ jLj 'iff ^ 

+ sign (0 — 0 + m) + sign (0 — 0 — m)\ , 


d0p(0) [-2|0-0| 


d 0 p( 0 ) [—2 sign (0 — 0 ) 


0 = zap"^) + + -M4>) 


4 r 2 

2p(0 + m) + 2p(0 — m ), 


(D.2) 


(D.3) 


(D.4) 


where equations ( |D.2 ),(D.3) and (D.4) corresponds to the first, second and third derivatives 
w.r.t. 0 . 

Now we are ready to consider solutions for different phases of the theory. In particular, 
we consider only solutions with no resonances and one pair of resonances. In the decom¬ 


pactihcation limit r —)■ oo we want to observe solutions (4.9) and (4.27). Unfortunately it is 


very hard to observe solution for the case of arbitrary number of resonance pairs n. However 
obtaining particular cases of n = 0 and n = 1 will strongly suggest that solutions we have 
found in section [| are right. 


D.l Solution with no resonances (n = 0) 

We start with the case for which we do not expect any resonances. This happens when 
2 a < m, where a is the position of the support endpoint as usually. In this case p(0± ?n) = 0 


as 0 ± m appears to be outside the support of the eigenvalue density. Hence (D.4) can be 
written in the form 


4/(0) - 4p( 0 ) = 0 . 

General solution to this differential equation is given by 

p(0) = C cosh(r0), 


(D.5) 


(D. 6 ) 


where C is an integration constant. Notice that we omit sinh(r0) due to the symmetry 

p{4>) = p{-4>)- 


Applying the normalization condition (3.14) to (D. 6 ) we obtain 

r 

Q — _ 

2 sinh(ra) 


(D.7) 
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Finally, substituting this solution into (D.2) we obtain the following condition for the position 
of the support endpoint a 


8tF , 

m - — I r — ar + coth(ra) = 0 , 


(D.8) 


To summarize, we have found that solutions for the eigenvalue density of A/" = 1 SYM with 
massive adjoint hypermultiplet and no resonances is given by 

r 


P(<t>) = 


2 sinh (raj 
= 0, 


cosh(r0), 


< a , 
> a, 


(D.9) 


and the eigenvalue support endpoint is defined by (D.8) 


Equation (D.8) can not be solved analytically. However in the decompactihcation limit 


r —> oo it simplifies a lot. Condition (D.8), defining the position of the support endpoint, 
reduces to 


m — 


87 ^ 

t 


(D.10) 


which exactly reproduces (4.10) found after neglecting the subleading terms from the very 


beginning. The eigenvalue density (D.9) becomes highly peaked around 0 = ±a with the 
width of the peak ~ In the decompactihcation limit r —> oo (D.9) can be considered 


as the sum of two 5-functions | (5(0 + a) + 5(0 — a)), because the peaks becomes infinitely 


high and narrow and integrate to one. So we see that (|D.9|) in the decompactihcation limit 
coincides with the solution (4.9) found in section [fj 


We can also compare the solution found here with the results for the numerical simulations 


of the full equation of motion (4.3). In particular in Fig,9(a) we show the results of this 


numerical solution with the orange dots. On the same plot we show the analytical solution 


(D.9) with the position of the cut endpoint found from (D.8) numerically. As we see analytical 


solution is consistent with the numerical results. 


D.2 Solution with one pair of resonances (n = 1) 

Now we can consider what happens with the increase of the coupling, when the support 
length becomes large enough to contain one pair of resonances, i.e. m < 2a < 2m. To 


solve equations (D.l) we should consider three intervals separated by the positions of the 
resonances 

Pa{4>)i m — a < 0 < a ; 
p(0) = p_b(0), — m + a < 0 < m — a; 

p C '(0), — a<0<— m + a; 


(D.ll) 


Symmetry requirement for the eigenvalue density p(0) = p(— 0) implies following conditions 
for the densities on the intervals A, B , C 


Pa( 0) = Pc(—0), Pb(0) = Pb(-0) • 


(D.12) 
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Figure 9: The eigenvalue density p(0) calculated with the following value of the parameters: 
r — 10, m — 5, , N — 200. The orange dots show results for the numerical solution, while 
the dashed blue lines show the solutions (D. 6 ) and (D.17) on (a) and ( b ) respectively. 


Now we consider (D.4) on the intervals A and B in order to find general solutions for pa and 


Pb respectively. Then we substitute these general solutions into (D.2) and (D.3) in order to 


fix integration constants together with the position of the support endpoint a. 

Interval A (in — a < 0 < a). Note that on this interval 0 + m > a so that p(0 + m ) = 0 
and also —a < 0 — m < a — m meaning 0 — m G C. Using these remarks and symmetry 


properties (D.12) we can rewrite (4.8) as the following delay differential equation 


ViM + jp'V™ -4>)~ 4r 2 P/i(^) + 2 r°-p A (m - <j>) = 0 


(D.13) 


General solution for this equation is given by 


Aa(0) = c 2 cosh (ar - y)) + ^3 sinh (^ r (V “ y)) > « 2 = ^ , ft 2 = ^ , (D.14) 


where coefficients a and 0 where found by substituting these solutions back into (D.13). 


This procedure is similar to writing out characteristic equation for ODE, however still a 
bit diff erent d ue to delays we have in our equation. In particular single exponents can not 
satisfy (D.14) and we should instead search for the solution in form of cosh (ar (0 — y)) and 
sinh (j3r (0 — y)). Shifts in the arguments can be seen as centering of the solution at the 
middle of the interval A which is made for the convenience and in principle can be adsorbed 
in constants and C 3 . 

Interval B (a — m < 0 < m — a). On this interval situation is simpler then on A, 
because 0 — m < —a and 0 + m > a so that both arguments are outside the support and 


hence p(0 ± m) = 0 leading to the following form of (D.4) 


4Pb( 0) -4r 2 ps(0) = 0, 


(D.15) 
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and its solution 


Pb(4>) = Ci cosh(r0), 


(D.16) 


where symmetry requirements (D.12) were taken into account. 

Interval C (—a < (J) < a — m ). Finally on this interval using symmetry properties (D.12) 
we obtain pc(4>) directly from the density pa(4>) given by (D.14). 

To summarize we have found 


0 - y)) + <C 3 sinh (fir (p - y)) , 


m 


a < 
a — m < 


< a , 

< rn — a , 


= C 2 cosh (yr 
p( 0 ) = Ci cosh(r 0 ), 

= C 2 cosh [ar (p> + —j j — C 3 sinh (j3r ^0 + — j j , — a < 0 < a — m , 

= 0 , 


(D.17) 


1 ,1 2 § 8 

101 > a, a = y ’ ^ = 5 ’ 


Now we should determine integration constants Ci, C 2 , C 3 and the position of the support 
endpoint a using normalization condition (3.14) together with equations (D.l)-(D.3), written 
out for different intervals. After some algebra we obtain the following system of algebraic 
equations 


Ci sinh(r(m — a)) + 2C 2 a 1 sinh (cur (a ——^ = - , 

Cs/r 1 cosh [pr (a - y )) + C 2 a~ l sinh (ar (a - y)) = ^ , 


4m 

IT 

4m 

"3 


2 a 

y 

2 a 

J 


—2 


4Ci cosh(r(a — m)) + 8C 3 0 2 sinh ( /3r ( a — — 


m 


167T 2 


t 


-h 2 r 


-2 


C 2 a 2 cosh (ar (a —— j 


+ C 3 0 


-2 


sinh (/3r (a —— j 


(D.18) 

(D.19) 

(D.20) 

167T 2 


t 

(D.21) 


Here in particular first equation (D.18) is obtained from the normalization condition (3.14). 
Second equation (D.19) is the result of the substitution of (D.17) into (D.3) with the as¬ 
sumption 0 G A. Finally, equations (D.20) and (D.21) are results of the substitution of 
(D.17) into (D.2) with 0 e B and 0 e A respectively. 

Equations (C.21)-(D.21) are transcendental and can not be solved analytically. Instead 
we can consider these equations and their solutions in the decompactihcation limit r —y 00 . 
Assuming that r(a — y) 1 and r(m — a) 1, or equivalently assuming that the radius of 
S 5 is large and we are far enough from the critical points, we can rewrite these equations in 
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the following form 


Cie r(m - a) + 2C' 2 cr 1 e Q ^ a -T 1 


= r, 


2 r 
~3~ 


) + C 2 a- 1 e“ r (°-^) 

2Cie r(m - a) +4C' 3 /3- 2 e /3r ( a -?) 

C 2 a- 2 e ar (“-?) + C 3 /rV r (“- 


4 in 

IT 

4m 


2 a 

T 

2 a 


-2 


167T 2 


-hr 

3 3 


167T 2 


(D.22) 

(D.23) 

(D.24) 

(D.25) 


Solving these equations we find the following expressions for the integration constants 
C x = re~ r{m - a) [1 - 2 (1 + /T 1 ) (a -1 + 3/T 1 + 4)] , 

C 2 = 2rae ar ( a 2 ) (l + (3 [a. 1 + 3/3 1 + 4) , ^ 2 g^ 


C 3 = 2 r^e" /3r ( a "^) 


^-(1+r 1 ) (a 


-1 


+ 3/T 1 + 4) 


The last unknown parameter of the solution is the support endpoint position a. To find it 
we substitute solutions for C 2 and C 3 into (D.25). Note that the last parenthesis containing 
C 2 and C 3 is of order r” 1 . Hence in the decompactihcation limit this term is irrelevant and 
we are left with the relation 

24vr 2 

a = 2m -— , (D.27) 

which reproduces the result (4.26) obtained in section[4j Moreover, we can reproduce solution 
(4.27) for the eigenvalue density from (D.17). To do this lets first notice that 5-function can 
be defined through the following limit 

5(0 —c) = lim f(r ), where f(r) = are ar ^~ c ' ) , 0<c (D.28) 

r—>-oo 

= 0, 0>c. (D.29) 

Indeed the function /(r) is peaked at 0 = c for large r and in the limit r —> oo this peak 
becomes infinitely narrow and infinitely high, while the integral of /(r) is normalized to 1 . 
Thus in the decompactihcation limit this function can be considered as the 5-function. We 
have already used this relation with a = 1 while considering decompactihcation limit of the 
solution (D. 6 ) with no resonances. 

The distribution (D.17) has similar exponential peaks at 0 = ±a and 0 = ±(m — a) and 
thus we expect to obtain the following eigenvalue density in the decompactihcation limit 

p(0) = c o 5(0 — a) + Ci<5(0 — m + a) + Ci5(0 + m — a) + c o 5(0 + a). (D.30) 


Coefficients Co and c 3 can be obtained from the solution (D.17) using the normalization 
prefactor in the definition of the 5-function (D.29) 

1 


Co = rfi (m-r 1 + -C 3 (fir)-' eM“-f) = - , 


1 


Cl = -C 2 (ar)- 1 e ar ( a ~f) - ±C 3 (/ir)” 1 ) + ^Cir" ± , 

2 2 2 b 


1 


(D-31) 
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where we have also used solutions (D.26) for the coefficients C \, C 2 , C 3 . As we see in the 
decompactification limit we completely reproduce solution obtained in section |4j 

Finally, we compare analytical solution (D.17) with the numerical results. We show this 


comparison in Fig, 9(b) The orange dots represent the results of the numerical solution, while 


dashed blue lines on the same picture show analytical solution (D.17) with the parameters 
Ci, C- 2 , C 3 and a evaluated numerically from the system of equations (D. 18)-(D.21). As 
we see from this picture analytical results coincide with the numerical ones perfectly. Small 
discontinuities at the resonances positions 0 = ±(m — a) are related to high sensitivity of 
the exponential factors in this solution as well as in the system of equations (D.18)-(D.21) 
which we solve numerically in the end. Hence these discontinuities are numerical artifacts 
and can be neglected. 
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